|gr-qc/0412114l 



Cosmological perturbations of brane-induced gravity and 
the vDVZ discontinuity on FLRW space-times 

in 
O 

o 

<N : Cedric Deffayet^ 

^ I lAP/GReCO, 98bis Boulevard Arago, 75014 Paris, France 

and 

Federation de recherche APC, Universite Paris VII, 
2 place Jussieu - 75251 Paris Cedex 05, France 



in 



(N 
> 



^ ' Abstract 

o 



We investigate the cosmological perturbations of the brane-induced (Dvali-Gabadadze- 
O • Porrati) model which exhibits a van Dam-Veltman-Zakharov (vDVZ) discontinuity when 

I , linearized over a Minkowski background. We show that the linear brane scalar cosmological 

^j^l perturbations over an arbitrary Friedmann-Lemaitre-Robertson- Walker (FLRW) space- 

time have a well defined limit when the radius of transition between 4D and 5D gravity 
is sent to infinity with respect to the background Hubble radius. This radius of transition 
r> \ plays for the brane-induced gravity model a role equivalent to the Compton wavelength of 

' the graviton in a Pauli-Fierz theory, as far as the vDVZ discontinuity is concerned. This 

well defined limit is shown to obey the linearized 4D Einstein's equations whenever the 
Hubble factor is non vanishing. This shows the disappearance of the vDVZ discontinuity 
for general FLRW background, and extends the previously know result for maximally- 
symmetric space-times of non vanishing curvature. Our reasoning is valid for matter 
with simple equation of state such as a scalar field, or a perfect fiuid with adiabatic 
perturbations, and involves to distinguish between space-times with a vanishing scalar 
curvature and space-times with a non vanishing one. We also discuss the validity of the 
linear perturbation theory, in particular for those FLRW space-times where the Ricci 
scalar is vanishing only on a set of zero measure. In those cases, we argue that the 
linear perturbation theory breaks down when the Ricci scalar vanishes (and the radius of 
transition is sent to infinity), in a way similar to what has been found to occur around 
sources on a Minkowski background. 
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1 Introduction 



Brane world models offer several interesting new ways to modify the gravitational in- 
teraction and mimic a 4D gravity theory from an intrinsically higher dimensional one. 
This can be achieved assuming that the (4 + n) dimensional bulk space-time has a spe- 
cial topology, with e.g. compact extra-dimensions [1], or geometry, being e.g. a patch of 
AdS^ space-time [2], and leads generically to modifications of the gravitational interac- 
tion at short distances. Variations of these models have been considered where gravity 
was modified at large (cosmological) distances, being mediated by a resonance of massive 
graviton propagating over an infinite volume bulk [3]. While the first models of this sort 
suffered from ghost-likc modes being related to the presence of a negative-tension dynam- 
ical brane (see [4] and e.g. the review [5]), the Dvali-Gabadadze-Porrati model [6], relying 
on an other construction, leads to similar properties for gravity, without suffering from the 
same drawbacks. This model has been applied to cosmology with the aim of explaining 
the observed late time acceleration of the universe by a large distance modification of 
gravity [7-9], with the exciting perspective that this can also be tested by solar system 
gravity tests [10,11]. Aside its interesting phenomenology, the Dvali-Gabadadze-Porrati 
(DGP in the following) model is also interesting from a more abstractly motivated view. 
It can serve as a consistent framework to study various long standing problems associ- 
ated with massive gravity. Here we are using this model as a tool to investigate the van 
Dam-Veltman-Zakharov (vDVZ) discontinuity [12] over a cosmological, but not necessar- 
ily maximally symmetric, Priedmann-Lemaitre-Robertson- Walker (FLRW) background. 
The essence of the vDVZ discontinuity is that the quadratic Pauli-Fierz [13] action for a 
massive graviton on a Minkowski background gives different physical predictions (like e.g. 
for light bending) from linearized general relativity whatever the smallness of the graviton 
mass. It is known, however, that this is no longer true on a maximally symmetric non flat 
background [14,15], when the graviton Compton wavelength (squared) is much larger than 
the background cosmological constant. The same is expected to be true for more general 
space-times with non vanishing curvature, where some curvature scale is supposed to play 
the role of the cosmological constant (square-rooted). Here we use some results obtained 
elsewhere [16, 17] for brane-world cosmological perturbations to investigate this issue in 
the DGP model, in which the graviton propagator has the same tensorial structure as in 
Pauli-Fierz theory. 

We show that the linearized DGP theory over an arbitrary FLRW space-time has a well 
defined limit when the radius of transition between 4D and 5D gravity is sent to infinity 
with respect to the background Hubble radius. This radius of transition plays for the DGP 
model a role equivalent to the Compton wavelength of the graviton in a Pauli-Fierz theory, 
as far as the vDVZ discontinuity is concerned. This well defined limit is shown to obey 
the linearized 4D Einstein's equations whenever the Hubble factor is non vanishing. This 
shows the disappearance of the vDVZ discontinuity for general FLRW background, and 
extends the above mentioned result for maximally-symmetric space-times. Our reasoning 
is valid for matter with simple equation of state such as a scalar field, or a perfect fluid 
with adiabatic perturbations, and involves to distinguish between brane space-times with 
a non vanishing scalar curvature and brane space-times with a vanishing one. In the last 
case, it is shown that in the limiting theory, perturbations are solely supported by the 
brane motion in an unperturbed Minkowski bulk. The question of the presence or absence 
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of the vDVZ discontinuity is thus first addressed here in the linearized theory. We also 
discuss the validity of the linear perturbation theory, in particular for FLRW space-times 
where the Ricci scalar is vanishing only on a set of zero measure. In those cases it is 
argued that the linear perturbation theory breaks down where the Ricci scalar vanishes 
(and the radius of transition is sent to infinity), in a way similar to what has been found 
to occur around sources on a Minkowski background. 

The paper is organized as follows. In the remaining of this introduction we remind 
aspects of the vDVZ discontinuity (subsection II. introduce the DGP model (subsec- 
tion ^21) as well as the solutions for the cosmological background that will be used later 
(subsection II. 3|) . In section [21 we recall some properties of cosmological perturbations in 
4D General Relativity (GR in the following). We then turn to discuss the cosmological 
perturbations of the DGP model (section • First, we introduce the general formalism 
that we will later use (subsection 13. Ij) , then we apply it to the study of the vDVZ dis- 
continuity as seen in the Newtonian potential exchanged between non relativistic sources 
(subsection l3.2|) . Next we turn to the case of time dependent perturbations, first by recall- 
ing the differential problem associated with this case f subsection 13.3(1 . then, by discussing 
FLRW backgrounds with non vanishing Ricci scalars (subsection 13.4(1 . before turning to 
time-dependent perturbations over a Minkowski space-time f subsection 13. 5(1 . We then ad- 
dress the case of a cosmological background with a vanishing scalar curvature, but a non 
vanishing Hubble factor (subsection 13.6(1 . before discussing the generality of some of our 
previous arguments (subsection 13.7(1 . In the following section (section ^ , we discuss the 
validity range of the linearized theory, in particular for those FLRW space-times where 
the Ricci scalar vanishes on a set of zero measure. Last, we summarize and conclude 
(section [3 • 

1.1 vDVZ discontinuity and Pauli-Fierz theory 

The van Dam-Veltman-Zakharov discontinuity [12] was first discussed in the framework 
of the Pauli-Fierz theory for a massive spin two field [13]. This theory can be defined by 
the following action 

5(^^)'"^) = 4'ii + 4? + ^Mc, (1) 

(2) 

where is the Einstein-Hilbert action 

^EH = -ttV / d^xV^ (r - 2A(4)) , (2) 

truncated to quadratic order into a small metric fluctuation h^i, , considered as a dynamical 
field over some background metric g^i, solution to Einstein's equations (with g^i, = g^y + 

h^y and k^^^ being the inverse squared reduced 4D Planck mass), S'pp'' Pauli-Fierz 
mass term, and Smc is the action defining the coupling between the graviton and matter. 
The Pauli-Fierz mass term is defined by 

4f^ I d'x^g h,yKp [r^r^ - r^f^) , (3) 
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where m is a mass parameter. On a flat background, this mass term gives a mass to the 
spin two field h^i^ in a ghost-free way. The couphng term is given by 

SMC = \j d'xV=I Vr''"' (4) 

where T^^ is the (conserved) matter energy momentum tensor. Note that by definition 
the action S'^^(-*)'°^ is just the quadratic action for a massless graviton coupled to matter 
and deduced from general relativity. Wishing to compare the massless (5(^(4)'°)) and the 
massive (5'*-'^(4)'™7^0)) theory, it is convenient to define the amplitude ^(^(4)'™) due to one 
particle exchange between two conserved energy momentum tensors Tfj_i, and S^i,. This 
amplitude is given formally by 

^(^W''"^ = ^jd''xS^"'{x)h^,{T){x) (5) 

= - nl^-^ j d^xd^x'S^'''{x)IL^,o.p{x,x')T''^{x'), (6) 

where h^^ (T) is the graviton field generated by T^j^, and is the graviton propagator. 

Let us first consider the case where the background metric g^i, is the Minkowski metric 
r]fj,i, and the cosmological constant A (4) vanishes. The dramatic observation made in refs. 
[12] is that, in this case, the two theories S^^'^^ and S^^''^^^^ give different physics whatever 
the smallness of the mass parameter m. Indeed, in momentum space, one has (j/^ being 
the 4-momentum) 

2 tt(™=0) ^ , 1 1 , /TX 

P = 2^t^aVi^P + ^ri^pVua - ^VnuTlaP + ■ ■ ■ (7) 

/ 2 , 2\ tt(™7^0) 1,1 1 , 

[p +m ) U'^^^p' = -ri^aVup + ^r]^i3Vua - ^V,^ur]ai3 + ■■■ (8) 



where momentum-dependent (and gauge-dependent in the massless case) terms have been 
omitted in the right hand side. The diff'erence in the coefficient of the third term in 
the right hand side of the above expressions is responsible for the discontinuity. Note 
indeed that in the massive case, this coefficient is mass independent. Considering e.g. non 
relativistic sources separated by a distances sufficiently small with respect to the graviton 
Compton wavelength, the amplitude due to the exchange of a massive graviton is given 
approximately by taking the m —>■ limit in ^(o>»"7^o)^ that we note ^(o>"*^o). One finds 

^(0,m-»0) ^ l^(0,m=0)^ ^9^ 

3 

so that the massive amplitude stays different from the massless one, whatever the smallness 
of the graviton mass. This translates into a similar discrepancy in the non relativistic 
potential between the two same sources, the potential of the massive theory being larger 
by a factor 4/3. This extra attraction can be attributed to the exchange of the helicity zero 
polarization of the massive graviton (which has 3 more polarizations than the massless 
one). This difference can be nullified by redefining the Newton constant of the massive 
theory with respect to the massless one, assuming, e.g., that one measures the Newton 
constant by some Cavendish experiment (indeed, if one does not do such a rescaling, 
the Newton constant of the massive theory would be given by 4/3 x k^^^/Stt = k^^^/6tt). 
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However, with such a rescahng, the discontinuity wiU then reappear in other observables, 
hke the hght bending. The latter will then be 25% smaller in the massive case than in the 
massless one, which is much too large to be compatible with current measurements of the 
light bending by the sun [12]. 

The situation is however very different when the background space-time is non flat. 
Indeed taking g to parametrize a de Sitter or anti de Sitter space-time with a cosmological 
constant A(4) , one finds that the previously defined amplitude A takes the form [14] 

^(A|.„m) ^ -„2^)|s'-'(Af +m2-2A,4))"'T„„ 

(2) 

where is the Lichnerowicz operator acting on spin-2 symmetric tensors [18]. Set- 
ting A(4^ to zero into this expression, one recovers the previous results for a Minkowski 
background. However, in contrast to the previous case, one has 

^(A(4)^0,m-+0) _ ^(A(4)7^0,m=0)_ ^-|^-|^^ 

So that the discontinuity disappears on a maximally symmetric background with a non 
vanishing cosmological constant [14, 15, 19]^. This shows that the limits m ^ and 
A(4) — do not commute; however the amplitude ^(^(4)'"*) goes smoothly toward ^(^(4)'°) 
when one lets m'^A'^^^ go to zero. 

As was alluded to here-above, the presence of the discontinuity can be attributed to 
the exchange of the scalar polarization (helicity 0) of the massive graviton. The latter 
is coupled to the trace of the energy momentum tensor of a conserved source. When 
the background is Minkowski, this coupling is mass independent and remains in the limit 
where the graviton mass is sent to zero. However, on a maximally symmetric space- 
time with non vanishing curvature, the helicity zero coupling becomes proportional to 
the graviton mass and disappears in the zero mass limit, allowing to recover the massless 
result. A similar disappearance can be expected to happen in more general cases where the 
background space-time has non vanishing curvature [20,21]. This is in particular supported 
by the study of Schwarzschild-like solutions where it has been argued some time ago by 
A.Vainshtein [22] that the vDVZ discontinuity could disappear non perturbatively, i.e. in 
the full exact solution of the theory. Indeed, as noted by Vainshtein, the Schwarzschild 
radius enters in combination with the graviton mass at the second non trivial order in 
the perturbation theory to open the possibility for a non perturbative recovery of the 
short distance behavior the standard Schwarzschild solution. In the case of a maximally 
symmetric background with non vanishing curvature, however, there is already a length 
scale (besides the graviton Compton length) at the first order of perturbation theory, the 
radius of curvature of the background space-time, and, as we saw above, this combines 
with the graviton mass to lead to a disappearance of the vDVZ discontinuity already at 
the linear level. 

One thus would wish to study the issue of the disappearance of the vDVZ discontinuity 
for general curved backgrounds. This is however plagued by possible inconsistencies. 

^However the theory is non unitary for A(4) > and < 2A(4)/3 [15]. 
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Indeed, the background metric g^u, discussed so far, was introduced as a completely 
extraneous field with respect to the "graviton" and one would need to have a fully 
non linear theory of gravity which linearized action would be given by (^. However 
there is no known non pathological such a theory. E.g. it has been shown [23] that 
starting from the full Einstein-Hilbert action, and going beyond the quadratic level in the 
expansion for the kinetic term, causes a ghost-like sixth degrees of freedom to propagate. 
Other non linear completion using two metrics (The so called strong gravity or bi-gravity 
theories [24,25]) have also been shown to suffer from the same drawbacks [25]. One could 
also think to start from a higher dimensional theory and truncate it in some clever ways 
to retain a finite number of massive graviton Kaluza-Klein states, such truncations are 
however also pathological [26]^. 

Here we would like to concentrate on a non linear theory which exhibit a vDVZ dis- 
continuity on flat backgrounds. This is the brane world model of Dvali-Gabadadze-Porrati 
(DGP) [6] that we will introduce in more details in the next subsection. So far, it has 
not been found to suffer from the same problems as the other non linear completions 
we mentioned, and in addition, it has known exact cosmological solutions [7]. The aim 
of this paper is to use those solutions as a consistent background to study the issue of 
the disappearance of the vDVZ discontinuity over cosmological space-times which are not 
maximally symmetric. 



1.2 DGP model 

The DGP model [6] we are considering is a 5D brane-world model with bulk gravitational 
action 

%) = ~:^72~ f d^^V9(^Ri5), (12) 

where i2(5) is the 5D Ricci scalar, and k^^^ is the inverse third power of the reduced 5D 
Planck mass. To account for the brane, one adds to this action a term taking care of 
brane-localized fields given by 

5(4) = J (fx^/g^^C, (13) 
where £ is a Lagrangian density given by 

In the above expression, is a Lagrangian for brane localized matter and R^^'^ is the 

Ricci scalar of the induced metric g^^} on the brane defined by 

4^= d.X^d^X^gfj,, (14) 

^Note also that the non perturbative disappearance a la Vainshtein of the vDVZ discontinuity has 
also been analyzed in some of the above mentioned non linear completion of action Q with a negative 
result [27,28] 
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X^{x^) are defining the brane position, where are bulk coordinates, and coordi- 
nates along the brane wor Id- volume^ . We also implicitly include in the action a suitable 
Gibbons-Hawking term [29] for the brane. It is always possible to choose a, so called 
Gaussian Normal (referred to as GN in the rest of this work), coordinate system where 
the line element can be put in the form 

ds'^ = dy^ + gjfjdx'^dx'', (15) 

and the brane is the hyper-surface defined by y = 0. In this coordinate system the 
gravitational equations of motion derived from the sum of actions (|12|) and l^l'A^ are simply 
given by 

Gfs = 4>Ay) fr(f ) - -^Gg)) (16) 



where is the 5D Einstein tensor, G^jD is the 4D Einstein tensor built from the induced 
metric (|14|) . and T^uf "^ is the energy-momentum tensor of the brane localized matter fields 
entering into C(^m)- 

Before discussing the DGP model defined by the equations of motion (|16)) . or the sum 
of actions (fT^ and (fT^ . let us first turn to a scalar toy model, proposed in the original 
reference [6], which captures some key features of the full DGP model. It can be defined 
by the action for the suitably normalized scalar field (j), 



[ d^xdy I -^dAc/^dU + ^^(y) J(4)<^ + -^6{y)d^^d''A . (17) 



1 . ^ 1 . . , l_ 

2 

(4) 

This model is simply one for a scalar field in a 5D flat bulk with a brane localized (in 
y = 0) kinetic term added for it, and a brane localized source term J(4). The equations of 
motion derived from this action read 

l^d^'dA + -^5{y)d^d>^\ ^ = ^<5(y) J(4). (18) 

In this model, the interaction potential between two static sources, separated by a distance 
r, interpolates between a 4D 1/r behavior at small distances and a 5D 1/r^ behavior at 
large distances, as shown in reference [6]. The crossover distance rc between the two 
regimes being is given by 

re = (19) 

^'^(4) 

The same behavior was found in ref. [6] to hold for the gravitational potential between 
two sources put on the brane in the theory defined by the sum of actions (|12|) and (|ll-{|). 

^Throughout this article, we will adopt the following convention for indices: upper case Latin letters 
A, B, ... will denote 5D indices: 0, 1, 2, 3, 5; Greek letters /i, v, ... will denote 4D indices parallel to the brane: 
0, 1, 2, 3; lower case Latin letters from the middle of the alphabet: will denote space-like 3D indices 

parallel to the brane: 1, 2, 3. 
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whenever the background bulk and brane space-times were both taken to be Minkowski 
space-times. However it was also found there that, from a 4D point of view, gravity is 
mediated by a continuum of massive Kaluza-Klein modes, with no normalizable zero mode 
entering into the spectrum. This being a consequence of the bulk being flat and infinite. 
As a consequence the tensorial structure of the graviton propagator was shown to be the 
one of a massive graviton and the model exhibits a vDVZ discontinuity. E.g. the Fourier 
transform of the amplitude, defined as in ©, between two conserved sources T^^, and S^u 
on the brane (which position is defined by y = 0), takes the form [6] 

h,AP,y = o)S'^-'{p) oc f , !r , (20) 

p-^ +rc p 

where inverse hats denotes Fourier-transformed quantities (and one has gone to the Eu- 
clidean space, so that p is simply the square root of the four momentum square). One 
notes above the factor 1/3 on the numerator similar to the one seen in the propagator 
(jHI, while the mass of the graviton appearing in (jH]) is here replaced by the "momentum 
dependent mass" r~^p, so that plays in the DGP model a role similar to the mass 
scale m of the Pauli-Fierz action ^. As discussed in the previous subsection, the vDVZ 
discontinuity would seem to rule out phenomenological applications of DGP model, if one 
trusts the linear approximation to compute observables such as the light bending from 
the sun. However, based on the fact that the exact cosmological solutions (found in [7] 
and recalled in the next subsection) do not exhibit any sign of the vDVZ discontinuity, it 
has been suggested that the vDVZ discontinuity of the DGP model was indeed an artifact 
of the perturbation theory [20] following the similar suggestion made by A. Vainshtein in 
the framework of Pauli-Fierz theory [22], and recalled in the previous subsection. This 
has been studied in different situations by expansions going beyond the linear order in 
the brane bending [27,30-33]. Those analysis, mostly concentrated on the case of spher- 
ically symmetric solutions on the brane, confirmed so far the original suggestion made 
in ref. [20]: for a spherically symmetric metric on the brane, the solution given by the 
linearized theory breaks down below a non perturbative distance, given by ^ 

rv = {r^crsf^^ . (21) 

This scale can be coined from the equivalent one found by A. Vainshtein to appear for 
a Pauli-Fierz model [22]. It depends on a curvature scale of the spherically symmetric 

^Indeed one can view as a corrective term in to the 4D action in the same way the Pauh- 
Fierz mass term © can be considered as a perturbation in of the massless action S^^jj. The vDVZ 
discontinuity is then stating that this perturbative expansion does not translates in a similar expansion in 
observables of interest. 

^Note however that all the spherically symmetric solutions studied so far are approximate solutions. 
This has several aspects. First, the recent work [33] suggests that the linearized theory is never a good 
approximation of the exact non perturbative solution of the brane (even if this work also finds the appear- 
ance of the scale 12111 as well as a non perturbative recovery of the massless tensorial structure below that 
scale). Second, the fact that the exact spherically symmetric solution is not known still allows a loophole 
to all analysis done so far. One could indeed argue that pathologies are arising in the exact solution, 
similar to those found in ref [28] in the context of bi-gravity. However the situation of the DGP model 
with respect to this issue is as good (or as bad) as in any other brane world models (like the Randall 
Sundrum model [2]) where the exact (stable) Schwarzschild-type solution on the brane, corresponding to 
the perturbative analysis, is not known. 
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metric, namely the 4D Schwarzschild radius rs- For distances much lower than it has 
been found that the spherically symmetric solution on the brane is close to a usual 4D 
Schwarzschild metric (with 4D, massless, "tensorial structure") so that there is no more 
discontinuity. In particular, one sees that diverges as Tc goes to infinity (this is similar 
to sending m to zero in Pauli Fierz action), the 4D parameters (n'^^-^, and the mass of 
the source) being fixed. This disappearance of the vDVZ discontinuity can be attributed 
to the fact that one mode entering into the gravitational exchange, and related to the 
brane bending (or the brane extrinsic curvature), has a cubic interaction which becomes 
important at distances lower than r„, as discussed in detail in [20]. This "strong coupling" 
is a blessing from the classical point of view, since it enables to recover the right tensor 
structure. A debate is currently going on on the quantum consequences of this strong 
coupling [34,35]. This will not matter for us, as far as the conclusions of this paper 
are concerned, since our discussion is here purely classical. Moreover we will mostly be 
concerned here with the linear perturbation theory. And, following the discussion of the 
previous subsection as well as what is suggested by the study of spherically symmetric 
solutions, ask whether already at the linear level, the vDVZ discontinuity disappears in 
the DGP model on cosmological background, when we let Tc go to infinity and keep the 
parameters of the 4:D part of the action, (|13() . fixed. We will come back on this "strong 
coupling" issue in section |3 

We thus need to look at the DGP equations of motion (|16() linearized over a cosmolog- 
ical background. That is to say we need to look at the so-called cosmological perturbations 
of the DGP model. We will only consider the case of scalar perturbations (as seen from an 
observer on the brane) since it is for those perturbations that the discontinuity is expected 
to show up, being related to the helicity-zero mode of the graviton. In the next subsection, 
we discuss some properties of the cosmological background we are interested in. 

1.3 Cosmological background for DGP model 

It is convenient to discuss the cosmological background in a so-called Gaussian Normal 
coordinate system with respect to the brane, where the background metric takes the form 

(is(5) = -n^{t, y)dt^ + a^{t, y)5ijdx'dx^ + dy"^, (22) 

where the 3D metric 5ij is a flat Euclidean metric (we will only consider in this paper the 
case of a spatially flat universe), and the brane sits at y = 0. One can further choose a 
time parametrization such that the function n is set to one on the brane (that we note 
n(fe) = 1, where the index (6) means here, and in the following, that the corresponding 
quantity is taken on the brane). In this case the induced metric (|14|) is simply given by 

rfs^4) = —dt'^ + a^^^^Sijdx^dx^ , (23) 

(where a^^) = a{t,y = 0)) and is of Friedmann-Lemaitre-Robertson- Walker (FLRW) form. 
Considering comoving observers to sit at fixed comoving coordinates on the brane, t is 
then simply the cosmological time on the brane. Accordingly with the symmetries of (|22j) 
the brane localized matter energy momentum tensor is taken of the form 

S'i = '5(y)diag [-p(M), P{M), P{M), PiM)) ■ 
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Solving the bulk Einstein's equations one gets then the expressions for n and a [7, 36] in 
the case of an expanding brane (i.e. when d(f,) is strictly positive): 

ait,y) = a(^h){'^-v\y\H}, (24) 
H + 

n{t,y) = l-r]\y\ — — — , (25) 

where H = a^p-^ /ci{b) is the brane Hubble parameter, and rj is the sign^ of the brane effective 
energy density p = P(^M) ~ '^H'^ / i^f^)- Note that the metric ()22() with a and n given by 
the above expressions, and y of a definite sign (i.e. considering only the y > or y < 
part of the space-time), is simply a reparametrization of 5D Minkowski space-time, with 
coordinate singularities where a or n are vanishing [7,37]. The time evolution of a(fe) is 
then obtained by solving the Friedmann equations of the DGP model, given by 

PiM) = -3H{P(M) + PiM)), (26) 
where we have defined the parameter T by 

r = ^. (28) 

The early time behavior of the cosmology defined by equations H26|l and ()27() is obtained 
by considering the limit T ^ 0. In this limit (which means that the Hubble radius 
is much smaller than the crossover radius rc), equation is identical to the usual first 
Friedmann equation and one recovers standard cosmology (with no sign of the vDVZ dis- 
continuity, as recalled above). At late time, however, one has deviations from standard 
cosmology, with the exciting possibility provided by the r] = —1 branch to yield an accel- 
erated expansion with vanishing cosmological constant [7-9]. Note that, for a fixed t] and 
a specified matter equation of state (or for given matter equations of motion), the Cauchy 
problem associated to the system ()26() and ()27() is well posed and the specification of the 
initial values of a(f,) and p(^M) lead to a unique solution for the scale factor on the brane. In 
particular, as will turn out to matter for the following discussion, one can choose a matter 
equation of state (albeit quite unconventional), and initial conditions, such that the brane 
background Ricci scalar -R(b) ' gi'^^n by -R^^) = 12H + 6H, vanishes at all times. 



2 Scalar cosmological perturbations in 4D general relativity 

Before turning to discuss the cosmological perturbations of the DGP model, let us here 
first recall some features of standard (4D) cosmological perturbations (see e.g. [38] for 
reviews) . We will do so in a way suited for comparison with the brane- world cosmological 
perturbations that we will discuss thereafter. We consider a given cosmological FLRW 
background space-time, solution of the Friedmann equations derived from 4D GR. We will 
work in the longitudinal gauge, where the (4D) linearized line element reads 

ds^ = -(1 + 2<^)dt^ + a^{t){l - 2^)5^jdx'dx\ (29) 

^77, which takes the values ±1, is equal to — e of reference [7] 
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where a{t) is the background scale factor, $ and ^ are the two left over scalar perturbations 
of the metric. For simplicity, we only consider a spatially flat cosmological background 
(so that the here-above 5ij denotes a 3D flat Euclidean metric). The linearized Einstein 
equations 

'^yT(^^=SG(^}, (30) 
then lead to the four scalar equations 



^UMm) = -2^-2-1>, (32) 

2 ^ — $ 

K-(4)^'^{M) = ^2 ' (33) 

kI^)6P(^m) = 2 (2- + 4 ) ^ + 2-i + 2§ + 6-^' + -^A($-^), (34) 

\ Qj Qj j Qi (X oQj 

where A is the Laplacian with respect to the comoving coordinates x*, the scalar compO' 
nei 
by 



nents 5p(^M)^^<l(M)^ ^'^{M)^ ^P{M) of the matter energy momentum tensor T^^!^^ are defined 



'^^0%/) = -^P{M), (35) 
STf^M) = di5q^M), (36) 

6T;^m) = ^^(M)5j+(A;.-i<5;.A)<5vr(M), (37) 

and are respectively the scalar energy density, momentum, anisotropic stress and pressure 
matter perturbations. In the above expression AJ- is defined by 

so that one has A = A^. Equations (|31II34() are respectively the tt, ti, traceless part 
of ij, and trace of ij components of the scalar part of Einstein's equations H3U|). The 
system ()31I34|) does not close in general. To close the system and solve for the evolution 
of cosmological perturbations once initial conditions are provided, one needs to add more 
equations, namely the matter equation of motion, or equation of state. A first simplifying 
hypothesis, that we will make in the following, is to assume a vanishing Jtt^j^,/), 

Stt^m) = 0. (38) 

This occurs in the case of a perfect fluid, but also when the only matter is a scalar field. 
Then equation H33|) leads to the equality between the two gravitational potential $ and 

^> = ^. (39) 

The potential ^ (or ^) can be determined by adding again an extra equation. E.g. in the 
case of adiabatic perturbations, the matter perturbations follow the equation of state 

SP{M) = CsSp(M), (40) 
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where c| is the sound velocity. Equation H40() can be used to ehminate 6Pi^m) ^^'^ ^P{M) 
from the Hnearized Einstein equation to yield the evolution equation for the gravitational 
potential <I> 



+ (4 + 34)^6 + 



2- + -,(l + 3c|) 



a 



r2 



0, 



(41) 



and all the other relevant quantities can be deduced from the knowledge of <I> as a function 
of time, as can be seen from equations (|3m34p . A similar reasoning holds for a scalar field 
which energy-momentum perturbations obeys the equation of state 



5P< 



(M) 



^P{M) 



+ 6H 



(42) 



where cp^^' is the background value of the scalar field. Equation 1)42(1 then yields in a similar 
way as for adiabatic perturbations, an ordinary differential equation for Note further 
that ()41[) is a second order differential equation so that one needs as initial conditions the 
values of $ and <I>. This will further specify the initial values of all the matter perturbations 
thanks to equations , , , (|in)) , and which can then be regarded as (initial) 
constraints. Last, we note that when one considers non relativistic matter sources on a 
Minkowski background (where a is constant with time and is set to one), one finds from 
equation (|^T|) the usual Poisson equation for the gravitational potential, reading with our 
normalization 



^(4) 



(43) 



3 Scalar cosmological perturbations in the DGP model 

3.1 General formalism and Mukohyama's master variable 

The most general scalar perturbations of the line element (|22|) read 

/ -n2(l + 2A) a^diB nA,. 

a2 [(1 + 2lZ)6ij + 2dfjE 



gAB 



■'I 
nA. 



y 



^y 

a?diBy 
1 + 2A 



(44) 



yy 



out of which one can define four gauge- invariant quantities. Following [40], we define them 
as A, Ay, Ayy, aud TZ given by 



A 




Ay + 



+ 



n 



n 



n 2- 
a, 



A 



yy 
TZ 



Ayy + {a'^ayY, 

- i_ aa _ 
IZ + aa ay ^a, 



(45) 

(46) 
(47) 
(48) 
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where a dot means a derivative with respect to t, a prime means a derivative with respect 
to y, and a and ay are given by 

a = -B + E, 

ay = —By + E' . 

In an analogous way, one can define a gauge invariant perturbed brane position ^ as [40] 

^ = ^-a^ay, (49) 

where ^ denotes the hnear perturbation of the brane position around the background value 
y = 0. The gauge invariant variables A, Ay, Ayy, R, as well as ^, are invariant under the 
scalar gauge transformations 

t t + 5t, (50) 
x' x' + 5'^dk6x, (51) 
y ^ y + 5y, (52) 

where 5t, 5x and 5y are arbitrary functions of t, x*, and y. The next step would be to 
introduce the decomposition (|44j) into Einstein's equations. However it will turn out more 
convenient for us to use rather Mukohyama's Master variable formalism [39], that we now 
recall. This can be done here, since the bulk background cosmological space-time of the 
DGP model is a slice of 5D Minkowski space-time. 

Mukohyama showed indeed that all the linearized scalar Einstein's equations over a 
maximally symmetric background bulk can be conveniently solved introducing a master 
variable Vt which obeys a PDE in the bulk, the master equation. The latter, when Vl has a 
non trivial dependence in the comoving coordinates x*, reads in the GN coordinate system 
(El 

0. (53) 



no? ) a? \ G? 

In the rest of this article, we will implicitly consider all the perturbations as Fourier 
transformed with respect to the a;*s, in order to do a mode by mode analysis. In particular 
H53|) can be rewritten as 

V/sSl = 0, (54) 

where Pa is a second order hyperbolic differential operator acting on y and t dependent 
functions (in the GN system), and A is understood to be replaced by —k'^, where k is the 
comoving momentum. Equation H53() is then only valid when k'^ does not vanish [39] , which 
is the only case of interest as far as cosmological perturbations are concerned. The gauge 
invariant scalar perturbations ()45I48() can be expressed in term of Q as follows [39,40] 

A = -J-(2n" + irn--^n---n'), (55) 

6a V n-^ n J 

Ay = —i^-'- -Q.-\ , (56) 
an \ n J 

Ayy = (VL" + 4^17- - - 2-n'\ , (57) 

6a V n J 

n = ^(n" + ^^- + (58) 

ba \ n-' n 
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Let us now turn to the brane. One can choose a gauge where the brane sits in y = 
and the perturbed induced metric on the brane has the longitudinal form (|29|) (with 
the brane scale factor a^f,) replacing a{t)). The brane localized matter perturbations are 

decomposed as in (jSSEZI), where T^if^ is now the energy momentum tensor of the brane 
localized matter fields. We will assume in the rest of this work that the matter anisotropic 
stress perturbation (Jvr^Af) vanishes. The Israel's junction conditions [41] together with the 
master equation H53() then enables one to express all the scalar degrees of freedom as seen 
by an observer on the brane, that is to say (^/0(a/), ^(1{m)^ ^P{m)-, ^ and as functions of 
the master variable O and its derivatives, as follows [16] 



^ = 4^{^^A(o,o)+^^^^-C(to) + ^^-C(l^^ (59) 

^ = ^^{^cl^o,o) + Hnxl,^ + nxl,^ + Hn%^,^^^^^, (60) 
2 H"^ j An p A^n p . p .. p An - p 

f^{A)^P(M) = I + Jp^^/\-^{Qfi) + '^(1,0) + ^ '^(2,0) + ;^^'^A{2,0) 

+^C(m + + w^^mi) + ^'^(i.dI^,) ' (61) 

+ (62) 

2 H"^ (An p p An- p p An p 




'jp^ilfl) + Hn'c^Q]^) + n 'c^ 1) + -^C(^ , (63) 

where the coefficients C are given in the appendix and our normalization is chosen so 
that they are dimensionless^ . We note here incidentally that despite the assumed vanishing 
of the matter anisotropic stress, equations H59|) and (|6U|) show that <I> and ^ are not equal 
in general. In fact, in the gauge considered here, the difference between <I> and ^ is related 
to the gauge invariant brane position by^ 

I = r,{^-^>). (64) 

Let us first apply this formalism to the case of perturbations over a Minkowski brane, 
and show how the vDVZ discontinuity does appear there in the Newtonian potential 
between non relativistic sources. 



3.2 vDVZ discontinuity on Minkowski background, as seen in the New- 
tonian potential 

Setting H to zero, a(ft) to one, in equations H59ll63p . and keeping Tc finite, one gets the 
following expressions for the gravitational potentials, the matter density, momentum, and 

*Note that the coefficients C axe given here in full generality for what concerns the sign of 77; while they 
are only given for the case r) = +1 in reference [16]. 

^This can be easily seen from the formulae given in ref. [16]. 
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pressure perturbations over a Minkowski brane as function of the master variable 



$ = -Af](,)--^,), (65) 



^> = ^A0(,), (66) 



iii)^P{M) = \[^^^- /^n - + , (67) 

^ I ) (b) 

4)^^(M) = -i|AO-0- + ij^-'j , (68) 

L J (6) 

K(4)^^(M) = l{Aa- -n - + ^n -'] . (69) 



2 



(b) 



As we noted above, $ is not in general equal to ^, this remaining true even if one sends 
Tc to infinity (after sending H to zero) as will be seen more explicitly below. It is the trace 
of the presence of the vDVZ discontinuity in the DGP model (and can be translated into 
the presence of a non vanishing anisotropic stress for the Weyl's fluid, which remains in 
the limit where rc goes to infinity on a 4D Minkowski background). It is interesting to 
recover from equations (|65II67|) the results obtained in [6] , and reminded in section 11.21 
for the Newtonian potential between two non relativistic sources. Consider indeed a non 
relativistic source which only non vanishing component of energy momentum tensor is 
(5/9(A/). Neglecting the time derivatives of the perturbations, equations (jH^ and read 

^ = ^^{b), (70) 

where we have defined ip hy Acp = Ail. The master equation ()53() together with equation 
((^7|l can be recast in the unique equation 

^ (Aif + c^") + -^6{y)A^ = ]-S{y)Sp^M) ■ (72) 

'^(5) ^^(4) ^ 

The above equation is the same as H18() and was solved in ref. [6]. Considering a localized 
source, then 'f[b)y ^ given by equation (|72l) . was found to vary for distances to the source 
much smaller than rc as a usual four dimensional potential verifying (at small distances) 
the standard Poisson equation 

^^Pib) = -Y^P(M), (73) 

while at larger distances it was found to vary as a 5D potential (see [6]). One then recovers 
from the above discussion the two-folded manifestation of the vDVZ discontinuity: keeping 
rc much larger that the distance to the source, one sees first from equations ()7U() and H73|) 
that the gravitational potential on the brane $ is renormalized by a factor 4/3 with respect 
to the usual one. This has been discussed below equation @. Secondly, one has ^ = 
which can also be derived from the amplitude (|20|) and be seen to match analysis of the 
Schwarzschild solution in the DGP model [6,20,27,31,32,42]. 
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3.3 General cosmological background: differential problem and well- 
posedness 

When dealing with a more general cosmological background, the problem we are interested 
in becomes the following: given specified initial conditions (that one has to provide both on 
the brane and in the bulk), say at a given cosmological time to, what are the perturbations 
as seen by a brane observer at a later cosmological time t > to? To answer this question 
for a given mode k requires solving the 2D PDE H54|) with initial (Cauchy or also possibly 
characteristic) data provided along some initial curve in the 2D {y,t) plane, and, given 
the hyperbolic nature of the problem, some boundary condition for 0, along the brane. In 
the following we will denote by C(7) the initial curve, and by C(fe) the curve representing 
the brane trajectory in the (y, t) plane. Those two curves are bounding a domain D where 
we need to solve the master equation (j^Sl- ^(7) and C(fe) intersect on the brane at some 
point O which represents in the (y, t) plane a space-like 3D initial hypersurface of cosmic 
time to- 

In general there is not a clear way to give a boundary condition on the brane, since 
this requires solving the equations of motion for the perturbations of brane localized fields, 
which in turn are intricated with the metric perturbations. However in the simplest case 
of matter with vanishing anisotropic stress (for which H59I63() hold), and simple equation 
of state, one can obtain easily such a boundary condition [16,43]. Indeed, for the case of 
adiabatic perturbations, substituting the expressions (|^ for 5P(^m)-, and (|HT|) for 5p(^M) in 
the equation of state (|ln|) yields a boundary condition on the brane of the form [16, 17,43] 



i^^rP ^HQCP ^^^'rP ^Q-rP ^^^"rP 

i -^'-A(0,0) + '-(1,0) + ;^^'-A(l,0) + '-(2,0) + ;^^'-A(2,0) 
+ -^<^(^,0) + 'Jp^{4,,0) + ^^'^{0,1) + ^ '^(1,1) + 'Jf^i2,l) 



2 I ^"" /^p _L '-^ /^p _L un nP 4- o rf 4- rP 
'^'-S S '^A(0,0) t72„4'"A2(0,0) ""(LO) " '"(2,0)"'" u2„2^/l 



^ '-(3,0) ^ -""'-(0,1) ^ ^^2'-A(0,l) ^ " (1,1) J 



(b) 

"^"-A(0,0) ^ jj2q4"-A2(0,0) ^ "-(1,0) ^ " "-(2,0) "^ ^2q2'-A(2,0) 

(74) 

The latter can be recast in the form 

= Y.^r,o)dl^ + Y.^r^)^l^', (75) 

r=0 r=0 

where the (3 are time dependent coefficients known from the background cosmological 
solution, and can be expressed as linear combinations of the coefficients C. A boundary 
condition of identical form is found in the case where the only matter on the brane is 
a scalar field, substituting the expressions for 5i-'(j\/), for 5q{^M)-, and (|HT|). into 
equation (|42|) . In the following we will restrict ourselves to the here-above mentioned cases 
where a boundary condition of the form H75|) holds. 

Equations (|53|) and (|75|) are then all what is needed to solve for the evolution of DGP 
brane world cosmological perturbations once initial conditions are supplied in the bulk. 
Those two equations play for the brane world cosmological perturbations of the model at 
hand, an equivalent role to the one played by equation H41j) for 4D adiabatic cosmological 
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perturbations of a perfect fluid. One might worry about the non standard form of the 
boundary condition (|75j) which involves derivatives of the master variable along the brane 
(such a boundary condition has been called non local [43]) . However, one can recast 
the differential problem defined by equations (|53|) and (|75|) in a standard form [17]. This 
is done as follows. First we introduce the characteristic coordinates X and Y of the 
differential operator T>^, such that the master equation is rewritten as 

dxdyn = ^ + J^ (76) 
= V{X)dYn + WA{X)^}, (77) 

where the coefficients V and W/\ are functions of X (and also of the wave number k, as 
far as Wa is concerned) defined by the first line One then defines the 15 unknowns 
,,), with 0<r<4, 0<s<4, and 0<r + s<4by 

u(r,s) = d'xdpn, (78) 
and recast the differential equation (|77j) in the following linear system 

dxU^r,s) = U{r+l,s), (79) 

for r + s < 4, 

dxu(^r,.s) = d^xdy^ (^^^(0,1) + ^An(o,o)) , (80) 

for r + s = 4 and s > 1, and 
^^^(4,0) = d\ (Fu(o,i) + TyAn(o,o)) • (81) 

The boundary condition (|75|) then takes the form 

^(4,0) = X! «(r,s)^(r,s)> (82) 
(r,s),r+s<4 

where a(4^o) vanishes by definition and the others a[r,s) ^-i'^ known functions of the param- 
eter t along the brane (the cosmic time) that can easily be obtained from the expression 
(|75|) . This rewriting of eq. (|75j) is possible as long as the coefficient of the highest deriva- 
tive of $7 appearing in ((7SJ), namely Vl "" for the cases considered here, does not vanish. 
The latter coefficient is given by 

2 + T 

2., (T + mjjL\ ^^^^ 



2a(h\ K 



It does not vanish for a sufficiently large class of cosmological backgrounds. The system 
H79ll8ip is a linear first order system written in its so called characteristic normal form, 
with a boundary condition (|82)) of standard form on a non characteristic curve C(b) (the 
brane is non characteristic, since it follows a time-like trajectory). Standard theorems on 

^'^Note that the differential operator 17611 is singular in X = 0, which corresponds to a coordinate 
singularity (see [16]) in the characteristic (X, F, x*) coordinatization of 5D Minkowski space-time. For the 
rest of the discussion one can (but does not necessarily have to) always consider a domain T) which does 
not intersect the X = line. 
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linear differential systems then insure that the differential problem associated with the 
system (|79I81|) and the boundary condition is well posed^^ provided one gives initial 
data of Cauchy type on the initial curve C(/), that is to say the values there of all the 
variables U(^r,s)j with r + s < 4. Under this condition, the problem has a unique solution 
in the domain D, and this solution depends continuously on any parameter which initial 
and boundary data may depend on. 

However, because obeys in the bulk an hyperbolic PDE of order two, one cannot 
freely specify on C(/) all the U(^r,s) an independent way without encountering inconsis- 
tencies. E.g., considering first the case where C(/) is non characteristic in the vicinity of O 
(and remains so all along the segment of the initial curve where we want to specify initial 
data), then it is enough to provide on the initial curve the value of 0, and of its normal 
derivative with respect to C^j^ in order to know there the initial values of all the Uj-^ with 
r + s < 4. This follows directly from the fact that the initial curve is a non characteristic 
curve for the PDE (|53() . On the other hand, if the initial curve is a characteristic in the 
vicinity of O, then, in order to know the initial values of all the (with r + s < 4), 

it is enough to specify on C(7) the value of as well as to specify in O, and only there, 
the values of the derivatives of normal to C(/) of order lower or equal to four (see [17] 
for more details). Note that once the values of the U(r^s) known on the initial curve 
C(/), and thus also at O, then all the initial values (i.e. at the initial time to) of the brane 
variables 5p(^M)-, ^Q{M)^ ^P(M) are known through equations H59II63|) . So that if one 

wishes to study a specific case where those initial values are specified, one should take 
care to choose the initial values of the (or rather, as we just saw, those of and its 
derivatives) accordingly. Moreover, as it will turn out to matter for the discussion in the 
next subsection, one should also demand that, in all case, the chosen initial values of the 
tt(^ 3) are compatible with the boundary condition ()82() . This is a standard requirement 
for a hyperbolic problem with a boundary. 

3.4 Disappearance of the vDVZ discontinuity on FLRW background 
with non vanishing scalar curvature 

We are now in a position to discuss the vDVZ discontinuity on a FLRW background. We 
want to study the behavior of cosmological perturbations of the DGP model where Tc is 
send to infinity with respect to a background cosmological length scale, that we will take 
to be the Hubble radius H~^. In this limit: i/ we want to see whether initial data for the 
brane quantities (5/9(jv/); ^Q(M)i ^P{M)-> can be specified arbitrarily close to the initial 

data for the corresponding quantities of usual 4D cosmological perturbations (or to say it 
in an other way, if the initial constraints of 4D GR can be fulfilled), and, ii/ we want to 
know if the evolution of those initial data differ, or not, from the one given by the theory 
of standard 4D cosmological perturbation as recalled in section [21 These two questions 
are the ones an observer living on the brane would like to answer in order to compare in a 
cosmological context the predictions of the linearized DGP model to the ones of linearized 
4D GR. 

Thus, the limit we will subsequently consider is Hrc — > oo, that is to say T ^ 0, 
keeping H, k^^^ , as well as the other 4D background quantities P(m) and -P{j\/) , finite. This 

^^Note in particular that there is one boundary condition 18211 for one in-going characteristic X = 
constant curve in the domain T) of interest, as it should (see e.g. [44]). 
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simply because we want to compare cosmological perturbations of the DGP model with 
the ones of standard GR in the same cosmological background. This limit will be denoted 
simply by T ^ 0. 

Let us first discuss the case of the background space-time, in a way appropriate for the 
discussion of cosmological perturbations. As we already explained below equation in 
the T — > limit, the brane Friedmann equations (|26II27|) go to the standard Friedmann 
equations 

P{M) = -^H{P^M) + P{M)), (84) 
3H^ = K(4)P(M), (85) 

obtained simply by making T = in the system (|26I27|) . For a specific equation of state 
for the background matter /9(a/) and P(m) i the systems (|84II85() and (|26I27() are both second 
order differential system in the scale factor a(b). This means in particular that one can 
freely choose the initial values of a(b) and d(b), at the initial time to- The brane Friedmann 
equations define subsequently a brane trajectory in the characteristic coordinates X and 
Y. The brane curve C(b) is indeed given by [16] 

X(b) = 0(6), (86) 
= (87) 

We denote by Cq,^^ any brane curve, which is given by the solution of the system H26I27() 
with a non vanishing T, and goes across some fixed point O in the {X, Y) plane, at which 
we fix accordingly the initial value for a^^) and d(fc) at the initial time to- We also define 
C('^)=° as the curve intersecting the previous ones in O, defined by the same initial values 
for 0(6) and d(fc) in O, but the system (|84I85|) . One sees that the limiting curve C'^f,-^^ 
is exactly C^f,)^^- This is only a geometric formulation of the fact that the early time 
DGP cosmology is identical to standard cosmology, as noted in [7], and that there is no 
discontinuity in the background cosmological solution. Note moreover that in the same 
limit as above, the background geometry stays perfectly regular (it is always a slice of 
5D Minkowski space-time) and the {y,t) coordinate system does not break down in the 
vicinity of the brane, as can be seen from equations ((211) and (|25j) . The conclusion of all 
this, is that in order to compare the outcome of cosmological perturbations in the DGP 
model and in standard 4D GR on a given FLRW background, it is enough to choose for 
the background equation of motion (|84I85|) and H26II27|) , the same initial values for a(fe) and 
d(b), and the same matter equation of state for both model. This is because the resulting 
trajectory C^^f^ can be made arbitrarily close to the purely 4D trajectory C^i,)^^ by taking 
a sufficiently small T. The difference between the two curves will then only influence 
the behavior of cosmological perturbations at the second non trivial order in T, which 
will be enough for our discussion. An alternative approach, which might be considered 
less physical since it requires an T-dependent equation of state, is to choose exactly the 
same curves C'^^f'^ and C^^^^ by taking a different equation of state for the background 
matter in the T 7^ case and the T = case. This leads to a slightly simpler discussion 
since the background brane space-time are identical in the T 7^ case and in the T = 
case. We will only consider in this subsection the flrst case: namely a given fixed equation 
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of state (or fixed matter equations of motion) for the background matter on the brane, 
with the understanding that this will result in an T-dependence of the domain T> over 
which we study the master equation. However we stress again that the two approaches 
are equivalent as far as the conclusions are concerned. We will moreover assume in this 
subsection that the limiting trajectory, Cj^^^, is such that the Ricci scalar, of the 

induced metric on the brane does not vanish on the cosmic time interval [io? ^max] over 
which we want to follow the evolution of the cosmological perturbations. This can always 
be done choosing initial conditions and the matter equation of state, and will also imply 
that the Ricci scalar of the background metric on the brane, for T 7^ 0, does not vanish 
on the time interval of interest, for small enough T. 

Let us now turn to the behavior of the perturbations themselves. The associated dif- 
ferential problem is entirely specified by the master equation . the boundary condition 
((75|) . and a domain V bounded by the brane C((,) and the initial curve C(7) on which we 
give initial data. Among those different elements, we have just discussed the behavior of 
C(f,) in the T ^ limit. As far as C^-^ is concerned, we pick a fixed initial curve in the 
bulk, independently of T (which we can always do). The curve C(/), which intersects all 
the C^^-j curves in O, is chosen so that it is non characteristic (a similar discussion with 
identical conclusions can easily be carried out when the initial curve is characteristic). 
Note that, when T does not vanish, the domain T) can be made arbitrarily close to the 
limiting domain bounded by the chosen initial curve and the limiting curve Cjj^^^. The 
master equation does not depend on T (see e.g. equation (|76|) ). so the only T-depending 
pieces that we have not discussed are the choice of initial data on the initial curve C(/) 
and the boundary condition (|75|) . Let us now turn to these last two issues. 

We first note one crucial result: assuming that does not vanish on [toj^max], it is 
easily seen, looking at the expressions of the coefficients C given in the appendix that 
all of those coefficients have a well behaved and bounded limit^^ when T ^ 0. This limit 
is continuous and the limiting coefficients, as function of H and its time derivatives, are 
simply given by the expressions of appendix 1^ taken in T = 0. A first consequence of this 
is that the boundary condition ()75|) has a well defined limit when T ^ 0. 

We turn then to the choice of initial conditions. For future convenience, we introduce 
first a new coordinate system r(i, y) and z{t, y) such that the curve C(/) is defined by a 
constant z, and that the r = constant curves are normal to C(7). For a fixed T, we can 
choose arbitrarily on C(7) the values of and of its normal derivative w.r.t. C(/). 
The only constraint these values have to obey comes from requiring their compatibility 
with the boundary condition (|75|) in O, which depends on T. Li addition, having in mind 
a comparison with cosmological perturbations of 4D GR, we wish to fix at the initial time 
Iq on the brane, that is to say at O, the values and (of ^ and <1> respectively). 
We thus impose three constraints on the initial conditions at O, which in turn impose 
constraints on the value of Vt, and of its derivatives at O, through equations (|75|) . (|59|) . 
and the derivative of H59() with respect to t. That these constraints can be simultaneously 
imposed at O, simply by choosing the values of 17 and on C(7) is easy to see. One way 
to proceed goes as follows. Using the fact the initial curve is assumed non characteristic, 

^^Provided one considers a background space-time where H and its time-derivatives are bounded on the 
time interval of interest. Note that one would have encountered divergences going as 1/T at any point 
where R^'^} does not vanish. This will be further discussed in subsection 13.61 
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one can express the derivatives Q", 0, " and 0, "" at O (that we denote respectively as , 



^to and O^^;-) as 



U — 2'^ ' 



0,, . 



where 0{n) denotes terms which are linear combinations of derivatives of ri, 5"r2 or d^d'^^l 
with r + s < n and s < 1. We then choose arbitrary values for 0, all along C(7). This 
means that all the derivatives of O of the form at known in O. In addition we pick 
some arbitrary value for at O. Then, inserting the expression (|88|1 in equation (|59|1 we 
determine the value at O of dzdr^ as a function of the chosen value for ^to- In ^ similar 
way we can determine the value at O of dzd'^Q from the chosen value for and the 
derivative w.r.t. t of equation (|59|) . Last the value of dzd^^l at O is determined from the 
boundary condition ((75|) at O. At this point the values at O of all the derivative dzd!^^, 
with r < 3 are known, and we complete the initial condition by picking any function dz^ 
along C(^j-j compatible with these values at O. Note that the above construction is possible 
because: i/ the highest derivatives of entering in the expression of ^ and ^, obtained 
from equation (|59() . and in the boundary condition H75j) are respectively fi", Q", and fi " ", 
and, ii/ the coefficients multiplying dzdr^, dzd^Q, dzd^Q in the above expression (|88|) . 

(p^ . do not vanish, nor diverge, because of the invertible character of the coordinate 
change z{t, y),T{t, y), and the non characteristic nature of C(7) (which insure that — z''^ 
does not vanish). More generally we call 17^^ and dz^Jj-j any set of initial condition on 
C(/) which are compatible with the three constraints mentioned above, for a fixed value of 
T. We will assume that the functions ^Jj^ and dz^Jj^ have continuous limits when we let 
T go to zero. We call S a set of initial conditions (for T G some interval ]0, Tmax]) having 
this property. Again inspection of the coefficients C and of equations (|HH|) . (|H!H) and (|!in|) 
show that such a set is not empty. To summarize, a set of initial conditions 

5 = ,T G]0,T^ax]} (91) 



has the property that its elements [il,Jjy dzi^Jj^j fulfill the boundary condition H75j) in O, 

give fixed values and for $ and <I> at O, and go to a well defined limiting initial 
condition when T is send to zero. We call the solution over T? to the master equation 
H53|) with the initial data (^QJj^ dz^Jj^^ chosen in a given set S and obeying the boundary 
condition ()75() . The here-above noticed fact that the boundary condition ()75() has a well 
behaved limit when T — > 0, as well as the fact that the differential problem we are looking 
at is well posed [17], shows that the solution 17 J has a well defined limit when T — > 0, we 
call this limit JlJ^*^. Once one knows this limit, one can compute from equations (|59|). 
CT . (IHTl) . (jnU and (jSni), the limiting value (as T ^ 0) at all time t G [to,imax] of the 
brane expressions "if, 5p(^M)i ^Q{M)j ^P(M)- The question of the vDVZ discontinuity 
then translates into asking how these limiting values, noted Spjj^^^, 6qjj^^^, 
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6P^j^ , compare with the corresponding quantities of standard 4D GR, than is to say the 

solutions of equations (|31I34|) with the same initial conditions and ^to- 

To answer this question, we first define new quantities, that we will denote respectively 
Sp^^j^jy '^^(Af)' ^^(M)^ taking the limiting values in T — > of the right hand side 
of equations ()59I63|) . E.g. is defined by 

^° = ^^|^CA(o,o) + ^^^C(to)+^^-C(to) + ^^^^^ , (92) 

where the index x=o means that all the coefficients C are taken in T = 0, but also that the 
quantities depending on the background brane trajectory (and thus also on T, as we saw 
above) are given by the limiting Cj^^^^ trajectory. This means e.g. that the H entering into 
H92|) is obtained by definition by solving the standard Friedmann equations H84II85() . Note 
further that the function entering into the above definition of '^^ , dp^j^^^ ^1{m)^ ^^'^ 
^^{M) ^° ^^^^ arbitrary. Now comes the remarkable fact: one can verify after some 
straightforward, but tedious, algebra that the expressions Sp^j^^-^, Sq^^-^, ^P^^m) 

verify identically the standard 4D perturbed Einstein's equations H31II34() with Sn^jif) = 0, 
for arbitrary Q and 0' (in particular equals '^^). This answers the question, since 
it shows that for any set if initial condition S (and thus any limiting solution fij^") 
the limiting expressions <I>^^°, 5pJ^^, ^iIm)^' ^-^(jif)^ exactly equal to the 

corresponding solutions of cosmological perturbations computed from 4D GR, with the 
same initial conditions. Thus the latter can not be distinguished from the former by a 
brane observer. In this sense, it shows that there is no vDVZ discontinuity on a FLRW 
space-time with non vanishing Ricci scalar curvature (with the restriction that our proof 
is only valid for matter with vanishing anisotropic stress and equation of state leading to 
a boundary condition of the form (f7S)) ). 

An other way to look at this result is to say that in the T ^ limit, the knowledge 
at O of <I> and ^ and of the brane matter equation of state are enough to compute on 
C(fc) the values of ^, ^, 5p(M)i ^Q(m)j ^P{m) at all time t S [toi^max]- So that the brane 
decouples from the bulk in this limit. And a Dirichlet boundary condition for the operator 
(f53|l is then obtained extracting from equations (|59|) . (|60|) . (|6T|) . (|62]) and l|63|) (and (|53|) ) 
the expression of 17 as a function of <I>, 5p(^M)-, ^Q{m)j ^P{m) (this is done in Ref. [16], 
the corresponding equation is obtained from equations (|109p and of this paper) . 

3.5 vDVZ discontinuity on Minkowski background, as seen in time- 
dependent perturbations 

It is instructive to contrast the case we have just discussed, with the case of time dependent 
perturbations over a Minkowski background^^. First, we consider the case of a constant 
H background (so that the background is either de Sitter or Minkowski). It is easy to 
obtain from appendix the expressions (|59l(i3|) corresponding to that case, by dropping 
all the time derivatives of H in the coefficients C. Doing that, the coefficients C are found 
to be polynomials of T and /x; e.g., for a constant H background, C^j-qq^ is given by 
/x(24 + 24T + 6T^). To obtain the expressions (|65I69|) one then sets H to zero in the 

^■^Time dependent aspects of the vDVZ discontinuity of the Pauh-Fierz theory have been studied in 
ref. [45]. 
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constant-i/ expressions. The essence of the vDVZ discontinuity, in the way discussed here, 
is that the hmits -ff — > 0, and Tc oo, do not commute. Indeed if one sets H to a constant 
value, and then let T go to zero, one obtains for example the following expression for 



1 



(An 



1 „ + 3Ha - a- + 27]Hn'\ , (93) 

which is equal to the expression of as it should be from our previous discussion, but 
stays different from (jbSj) and (|66() in the — > limit^^. 

We would like now to ask the same question as in subsection l3.4( but for time dependent 
perturbations over a Minkowski brane. First we recall that if we set (57r(jvf) and 5P(^m) to 
zero in equations (|31I34|) (with d = 0), or equivalently cs to zero in equation (|41j) . one 
finds that <I> equals ^ at all time and is given by a affine function of time. If one further 
demands that the brane is empty at the initial time to (that is to say that (5/0(Af) and 
vanish at initial time), then the initial values of $ and $ have to vanish thanks to the 
initial constraints H31|) and H32|) . meaning further that vanish at all later time. As is 
well known, there are no scalar gravitational waves in 4D GR. Let us now discuss the case 
of the DGP model, with a Minkowski brane, the same initial conditions, and an empty 
brane. The master equation (|53|) takes the form of a usual 5D d'Alembertian equation 

n - - n" + Pn = 0, (94) 

a solution of which can be taken to be simply the "zero mode" 

n = Asm{u;{t -to) +uJo) , (95) 

with a;^ = k^. The 5p(^M) defined from equation (|H7|) and this solution vanishes identically, 
as well as the (5g(M) ™d , so that this solution verifies the boundary condition for an 

empty brane. On the other hand the initial values and can be chosen arbitrarily 
by choosing the amplitude A and the phase ujq to verify 



tan Wo = (96) 

as can be seen from equation (|65j) . There are thus propagating scalar gravitational waves 
on the brane in the linearized DGP model. Those waves verify ^' = — as can be seen 
from equation (|66|) and (|95|) , in contrast with 4D GR, where the vanishing of 5'K(^m) implies 
the equality between <I> and ^ . If one wants then to stick with the reasoning of section 
13. 4| one can choose initial data on some fixed initial curve C(7) such that the solution 
H95|) with the choices H96|) will be the unique solution corresponding to an empty brane, 
initial values and , and those initial data. The initial data, so defined on C^j^ , will 
obviously be independent of rc, as is the solution (j^Sl), so that both will be unchanged 
in the rc — > oo limit. This provides an example that the recovery of solutions of 4D GR 
cosmological perturbation theory shown to occur on a FLRW brane with non vanishing 
scalar curvature, does not happen for a Minkowski brane^^. 

^■^Equation (1931 also illustrates the fact that, on a cosmological background, the perturbation theory 
depends explicitly on the phase 77 = ±1 of the background, this even for distance scales unrelated to rc, as 
first noticed in reference [10]. 

^®It is quite simple to find other solutions showing the same, with rc dependent initial data. 
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3.6 Case of a brane with non zero Hubble factor, but vanishing scalar 
curvature 



We now turn to the case of a background with a vanishing Ricci scalar Rjly Namely, 

we will assume that R^^-^ vanishes at all times in the time interval [to,tmax]- This can be 
achieved as discussed in subsection ll.31 In this case, the reasoning done in section ITU is no 
longer valid because the coefficients C all diverge as 1/T, when T is sent to zero. However, 
one can easily get rid of this divergence by considering the new variable Q defined by 

n = rn. (97) 

O obeys a PDE in the bulk that can be obtained from (|53j) and reads 

^ n ) +2h( -X.) n + 4H( -Xr) n- - ^aq - ( ^n')' = 0, (98) 



where H means a(fe)/a(b), as in the preceding. To obtain this expression, we have used the 
vanishing of i?^^ which implies that 

t = 

= 2HT. (99) 

It might seem that the above PDE depends implicitly on T through H. This is however 
not the case, since H is fully determined by solving the equation 2H^ + H = 0, once 
the initial values of a^^) and d(ft) have been provided. This in turn also determines the 
brane trajectory in the characteristic coordinates X and Y of the differential operator 
associated with equation (jUSJ. Note that the latter coordinates are the same as those 
previously introduced, since they only depend on the derivatives of second order in the 
differential operator considered. Substituting O by Tfi in equations ()59II63() one obtains 
expressions for <I>, Sqt^M)-, ^-^d 5P[m)-, as functions of Q. and its derivatives, with 

new coefficients C, instead of C. One can verify that those new coefficients have all now a 
well behaved limit when T is let to zero. In addition, one can show as in section EU that 
the limiting expression of \I', 5p(^M)i ^Q{M)i a-nd 5P(^m)^ defined as <I>'^, Sp^^j^^y Sq^^jy 
and 5P^]^^ , with C and Cl replacing respectively C and 0, verify identically the linearized 4D 

GR equations (|31II34|) . whatever the function Q. This means that the rest of the discussion 
follows as in 13.41 

The case discussed here is however quite different from that discussed in subsection 
13.41 Indeed, the limiting values of O and its derivatives are bounded in the vicinity of the 
brane, which means that those of and its derivatives vanish by virtue of the definition 
1)9 7|) . It implies that the limiting values of the four 5D scalar gauge invariant variables 
H451148() vanish; this because those gauge invariant variables are all expressed as linear 
combinations of 0,, its derivatives, and the background metric [39,43] (see equations (|f)5t - 
158(1 ). On the other hand the limiting value of the gauge invariant brane position ^ is finite, 
as can be seen from equation (|(i4|) . Thus, the limiting regime is describing a brane moving 
in an unperturbed bulk, that is to say here, a slice of 5D Minkowski space-time. So, the 
limiting perturbations of the induced metric on the brane gj^J , defined by equation (|14() , 
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should be solely supported by perturbations of the embedding functions X {x^^), with 
of equation ()14|) given by the background cosmological metric ()22p. One may think that 
such a situation can only occur for specific cases (see e.g. [46]). As a consistency check, 
one should be able to verify that one can indeed embed any scalar perturbation of the 
brane induced metric of the limiting form (that is to say obeying linearized 4D Einstein's 
equations) in a unperturbed 5D Minkowski bulk when the background brane Ricci scalar 
vanishes. Let us now show that this is the case. 

We start with a brane defined as a hypersurface of coordinates X'^(x^) in the space- 
time parametrized by the metric 1)221) . The background embedding function X^^{x^) are 
given by 

^(B)(2;^) = '^^^^ (100) 

so that one can write at linear order the X^(x'^) as 

X° = t + 6t, (101) 
X' = x' + 5'^dk5x, (102) 
= 5y, (103) 

where we have considered only embedding functions generating scalar perturbations on 
the brane world-volume. Defining then the linearization of the brane induced metric as 



-(1 + 2^W) ^l^d^B^^^ 



al,^d,B(^^ af^ (l + 27^(4))5,,+252.i^(4) 



(104) 



equation 1)141) leads to the following expressions for ^(4), 5(4), and E^^"! as functions 
of (5t, 5x and 5y 



h ^ n' 



^(4) = st + -6t + -6y, 
n n 



7^(4) = '±5t + -5y, 
a a 

5(4) = -t,6t + 6x, 



0? 



£;(4) = fe, 



where the functions a and n appearing in those expressions are given by those of equations 
H24II25() evaluated in y = 0. If one then insists that the induced metric has the form H29|) . 
one sees that one should choose bx and bt to vanish, and by to obey the two equations 

n' 

$ = —by, (105) 
n 

^ = --by. (106) 
a 

We wish in addition to impose $ equal to ^, but ^ otherwise arbitrary. We see from the 
system (|105I106() that this is possible only if the coefficient multiplying by in 



$ - ^» = 




(107) 
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vanishes identically . That is to say if 

/6 = H + 2H^ = 0, (108) 

where, to obtain this expression, we have used equations H24|) and (|25|) . Thus we conclude 
that arbitrary scalar perturbations of the brane can be supported only by fluctuations of 
the brane position if and only if the background Ricci scalar vanishes. This is then fully 
consistent with the calculations done in the first part of this subsection. 



3.7 Initial conditions and the limiting theory 

In the discussion carried in the previous subsections 13.41 and 13.61 we considered sets of 
initial conditions on C(/) chosen in a subset of all possible initial conditions, namely sets 
of the kind l|9H) . Those initial conditions were chosen in particular so that they lead to a 
well behaved master variable ^g^^ on the brane. One could however worry that choosing 
a less restricted class of initial conditions could result in a limiting theory different from 
4D linearized Einstein gravity. We would like to argue here that this is not the case unless 
the limiting Q (or some of its derivatives) is diverging, leading to a subsequent divergence 
of bulk metric pertubations. Indeed, one can show that 17 fulfills the following equalities 
on the brane [16] 

A2j1(^) = -6a^ {3H (^"^ + H^^ - a-^A^ + r]T-^Kl^)6p}^^y (109) 
A^^) = -a^H-^ ^^a-'^^H^ An- 2rjT-^Kf^-^{6p- 3H6q- A6tt)}^^^, (110) 

where we have defined 5p, 6q, 5tt and (for future use) 6P as 

Sp = 5p(^M) + ^U-h + -'A-2^^] , (111) 
Sq = '5g(M) + ^|-«J> + ^j , (112) 



_ $1 

I ''(4)" J (6) 

6P = (5P(M)-^<i (4- + 24)^ + 2-6 + 2§ + 6-^' + -^($-^')l .(114) 




On the other hand, 0, is expressed in terms of the gauge invariant scalar perturbations of 
the 5D metric, A, Ay, Ayy (defined in equations (|45ll47j) ) as [16] 

2aa' 
3^ 



A^r^ = aa' (4:A{A + Ayy) + Gaff^^iA + 2Ayy) - 3a^i,^a' Ay 
+6od^ft) (2d(^h)^y + a^' + 2a'Ji') 

+3ad(b)a' (2ii^b) (^ - Ay) + + a'Ay^ , (115) 



^^Note further that whenever H is non zero the same is true for a' /a, so that a non vanishing Sy can 
indeed generate a non vanishing "if perturbation, in contrast to the case of a Minkowski brane. 
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which is holding everywhere in the bulk (and thus also on the brane). Notice then that the 
right hand side of equations (|111II114'|) are the linearized 4D Einstein equations, so that 
the non vanishing of 5q, Sir and 5P measures the failure of the perturbation theory 
to match with usual 4D linear cosmological perturbations theory. Let us now consider, 
e.g. equation H1U9() in the limit T ^ 0. With a choice of initial conditions ()91() and a 
background with a non vanishing scalar curvature, we saw previously that 6p vanishes as 
T, when we let T go to zero. It is then easy to see from equation H1U9() that the only 
possibility for 5p to have a non vanishing limit as T goes to zero, is if A^il^^) diverges 
simultaneously^'^. Now if A^Q diverges on the brane, equation (|115() shows that some of 
the gauge invariant variables (or some of their first derivatives) have to diverge as well, 
leading to a breakdown of the linearized perturbation theory. A similar reasoning can be 
made starting from equation (|110() . or in the case of a background with vanishing Ricci 
scalar. 



4 Breakdown of the linearized theory 

So far we have only discussed the linearized DGP model. Namely, we have rephrased the 
question of the presence or absence of the vDVZ discontinuity over cosmological back- 
ground into that of knowing what is the limiting theory, as rcH ^ oo of linear DGP 
cosmological perturbations theory, when such a limit exists. This is a question one can 
ask whatever the validity range of the linear theory is. We would like now to question 
this validity range. As we will see, this can shed also an interesting light on the vDVZ 
discontinuity, and in particular on the difference between the cases with vanishing and non- 
vanishing background Ricci scalar. For reasons that will become clear in the following, 
and motivated by a similar discussion held for the metric around spherically symmetric 
bodies [32,35], we will concentrate on the brane extrinsic curvature X^t^ The latter is 
expressed in terms of the brane effective energy momentum tensor S^i, by the Israel's 
junction conditions [41], derived from equation 1)16(1 




(116) 



where the brane effective energy-momentum tensor is given by 

Note that this tensor vanishes identically if the 4D Einstein equations are obeyed on the 
brane. Inserting this equation in the trace of equation (|116() gives 

Ki'-> = Si (t,„, + . (118) 

where K = K^ji_yg^^y ^(M) = ^lif^^^y and ii(5) = —G^^^g^^-^ is the Ricci scalar of the 
induced metric. Moreover, a Gauss decomposition of the 5D Einstein tensor leads to the 

^^We remind that in this hmit, we always ask the 4D "observables" on the brane, such as $ or $ to be 
bounded. 



27 



equation (see e.g. [36]) 

R('^ = kC^^K^,^ - Kf)K^^. (119) 

Extracting the expression of Rq,) from equation (|118() and inserting it into equation H119|) , 
one finds the equation 

Ik(^) - + K^}K^^^ = kI.^T^m), (120) 

which relates the brane extrinsic curvature to the brane matter sources. Note that this last 
equation is valid in full generality. Let us now consider truncating the theory to linearized 
perturbations. In this case one can decompose the extrinsic curvature tensor as 

K^, = K\^)+kI^\ (121) 
where K^!^'' is the background value of ivT^jy , and K^i^ , the piece of the extrinsic curvature 



linear in the metric perturbations. Linearizing equation p20j) over the background metric, 
one obtains 

^<) + ^W = 4)^S' (122) 

where -f^(L) is defined by 

and is linear into the metric perturbations. Notice that if one goes beyond the 

linearized theory, the left hand side of equation (|12(H) contains necessarily the term 

Let us first discuss the case of a Minkowski background. In this case -RT^^ is vanishing 
so that one expects the linear perturbation theory to break down at least when the terms 
(|124|) . quadratic in the linearized extrinsic curvature, become of the same order as the first 
term on the left hand side of equation H122|l . This is what is found in references [32,35] 
for the case of a non relativistic point like source. The breaking down of the linear 
perturbation theory can thus be estimated to happen when 

or 

Kl'jlfc ~ 1. (126) 

In the case of a non relativistic point like source, we note that equations ()116() and H117|) 
show that xj^i^^ is of the order^® of TcG^S^^^. We should then have for dimensional reasons 



^^One could have thought that S^v is vanishing as r^T^, when one lets rc go to infinity. This is however not 
so. Indeed, because there is a vDVZ discontinuity on a Minkowski background, the 4D Einstein equations 
are not fulfilled in the limit of large Vc, and there remains a non zero contribution in S^^. 
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^(b) ~ ''^cfsr ^, where rs is the (usual) Schwarzschild radius of the source, and r is the 
distance to the source. Thus, from equation H126|) one concludes that the linearized theory 
is no longer valid for distances to the source smaller than r„, r„ being given by equation 
(|21|) . This is what is found in ref. [20,27,30-32,35] in agreement with the similar behavior 
found in the Pauli-Fierz theory [22]. Notice that r„ goes to infinity as Tc goes to infinity. 
This indicates in particular that, for a fixed source, and at a given distance to this source, 
the linear perturbation theory is never a good approximation as Tc tends to infinity. 

Let us now discuss the case of a cosmological background, and concentrate on the kind 
of non linearities 1)124(1 that are responsible for the appearance of the "Vainshtein" radius 

as recalled above. In this case, K^iu is no longer vanishing. Rather, from the identity 

Kl^J = {ldy9l!}}^_^, (127) 

valid in a Gaussian Normal gauge, where the 5D metric takes the form (|15() and the brane 
sits in y = 0, as well as from equations (|24I25|1 . one finds the non vanishing components 
of Kjju to be given by 



= V^^, (128) 
^^^^^ = -vHal,^6,,, 



4^)(^) = -vHal,6ij, (129) 



SO that is of the order H. This means in particular that K(^i^ is likely to dominate 

over the first term in the left hand side of equation ()122() . when one considers the limit 
T ^ 0. This can indeed be checked as follows. We want first to estimate the components 
of in term of 17 and its derivatives, as well as their dependence in re- To do so, we 
work the GN gauge (|15j) . where the 5D perturbed metric (|44|) is given by 



/ -n2(l + 2^) a'^d^B 



9AB 



a?d,B a2 [(1 + 2n)5^j + 2dijE | , (130) 

V 



and we have kept only scalar metric perturbations. A, B, TZ, and E^^. The following gauge 
conditions can be further imposed without loss of generality 

Eib) = 0, (131) 
Bib) = 0, (132) 

and the induced metric on the brane is given by the line element (|29|) with $ and ^ given 
by 

Aib) = (133) 
7^(b) = (134) 

(135) 



So that A, B, TZ, and E are the values of A, B, TZ, and E in the gauge considered 
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Israel's junction conditions H116() translates into the following set of equalities [16] 



^2 

= -f^(3'5P + 2M, (136) 
7^'(,) = ^4)(A57r-M, (137) 

2 

Using those formulae, as well as equations (|59I63() . it is easy to extract the Vc dependence 
of the linearized extrinsic curvature K^^^^-^ (as well as for K^J^^^ and Kj^^^^^). Let us first 
assume, as previously that the background Ricci scalar does not vanish on some time 
interval [to,tmax]- One finds then that, at the leading order as T goes to zero, -f^l^") is 
given by a sum of terms of the form 

a^,,,,r)dldlAPn, (140) 

where p G {0,1,2}, q G {0,1}, and r £ {0,1,2,3,4}, and Qi(p,g,r) are time dependent 
coefficients depending only on the background metric through the scale factor a and its 
time derivatives. Notice in particular that the linearity of K^i, in Tc (through k^^^), implied 
by equation ()116|) (or equations H136II139() ). has disappeared. This is because the right 
hand side of equation (|117|) is found also to depend linearly on at leading order, when 
considering a cosmological background; or, to say it in an other way, as we saw in section 
13. 4( because the right hand sides of equations (|1 11111131) vanish as T, as we let T go to zero 
and keep 0, (and its derivatives) bounded. This contrasts with the case of a flat Minkowski 
background, as seen above. Notice also that equations (|61II63|) lead to a form similar to 
that of Kjj^-^ for Tj-^^ . We know in addition that the limiting expression of T^^^-^ , as T goes 
to zero, can be made non vanishing, this because it is obtained solving the 4D linearized 
Einstein equations (|31I34|) . While in the same limit, we conclude from the form (|14()|) that 
K^j^^/vc is in fact vanishing, so that the dominant term in the left hand side of equation 

(|122|) is of order K^^^H. The next step is then to ask whenever the nonlinear term 
(|124|) become of the same order as the dominant linear ones. This happens now when 

Kfj>^H-' ~ 1, (141) 

which does not contain re- This means that when T goes to zero, one should find a 
modified, rc-independent, Vainshtein scale'^''. This is indeed what has been found in 
references [10,47]. The same reasoning and conclusion would be true in the case of a 
background with a vanishing Ricci scalar on the time interval \tQ,tmax\- In this case, the 
preceding discussion is unchanged once one replaces Q by $7. 

Last, we would like to ask what is happening if one considers a background with a 
Ricci scalar that vanishes only on some set of zero measure. Consider, e.g., a background 



I Thank Arthur Lue for discussing this point. 
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such that its Ricci scalar does not vanish on some interval [to, tmax[, but vanish at the end 
of it, in tmax- Then, fix some initial data in ^-nd in the bulk, following the procedure 
explained in section [3.41 and consider the limit Hrc ^ oo. If one does so, one will find 
that the master variable Q, as well as its derivatives, stay perfectly well behaved in the 
bulk^^. This is then also true for the bulk gauge invariant variables A, Ay, Ayy and 
TZ (see equations (|45I48|) ). and so also for the bulk geometry. On the other hand, the 
4D observable quantities on the brane, ^,"^,6p(^m)t^Q{m)^ and ^^(a/) have in general a 
bounded (and continuous) limit for t < tmax, but diverge in tmax- This is because of the 
divergence, underlined in section (|3.6|) . of the coefficients C, as T and i?^^ go to zero. 

The same derivation that lead to ()140() shows that K^j'^^^^ is also diverging in t = tmax- 
One thus concludes that the linearized theory is breaking down when one approaches this 
moment. A similar conclusion is reached if one considers the gauge invariant brane location 
^. As can be seen from equations and (|59l60j) . the limiting ^ is finite for t < tmax, 
but diverges at tmax- Those divergences constitute somehow a reappearance of the vDVZ 
discontinuity when one goes continuously from a background with non vanishing Ricci 
scalar to a background with a vanishing one. They are of course likely to be absent in the 
full non linear theory. 



5 Discussion and Conclusion 

In this paper, we discussed the theory of linear cosmological perturbations of the Dvali- 
Gabadadze-Porrati model. We distinguished the case of a brane background space-time 
with a non vanishing Ricci scalar from that of a brane background space-time with a 
vanishing Ricci scalar. In the first case, we showed that if one fixes at some initial time 
along the brane world-volume, the initial values of the gravitational potential $ and its 
time derivative as can be done in 4D general relativity, and if one further insists 
that the initial data to be provided in the bulk have a well-behaved (continuous and 
bounded) limit when we let the radius Tc of transition between the 4D and 5D regime 
go to infinity with respect to the Hubble radius of the background space-time, then the 
linearized theory on the brane does have a limit, and this limit is given by the theory 
of cosmological perturbations derived from standard 4D GR. This holds for the simplest 
case of matter living on the brane, namely a perfect fluid with adiabatic perturbation, 
or a scalar field. This result can be rephrased in saying that the boundary condition 
on the brane, derived from the matter equation of state, and that is usually "non local" 
in time [43], degenerates, in the limit considered, into a standard Dirichlet boundary 
condition that is obtained solving the linearized 4D Einstein's equations. So, in fact, in 
this limit, the brane decouples from the bulk. It contrasts with the case of a Minkowski 
background, and can be interpreted as showing, in a "Hamiltonian" sense, that the vDVZ 
discontinuity does disappear on a FLRW background in a similar way as it does on a 
maximally symmetric background with a non vanishing curvature. What we mean here 
by that, is that an observer living on the brane, and considering only scalar perturbations, 
thus, say, fixing some initial data as he would do in standard 4D cosmology, will not be 
able to distinguish the time evolution of these data in the DGP model with large enough 
VcH, from what he would have observed in standard cosmology with the same initial data. 

^^disregarding some possible singularities arising due the hyperbolic nature of the problem. 
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The same result is found to hold in the case of a brane background space-time with a 
vanishing Ricci scalar, with however the important difference that the initial data in the 
bulk are to be such that their limiting values vanish. This contrast with the previous case 
where the only constraint those bulk initial data have to fulfill are those implied by the 
choice of initial data on the brane. Thus, in the case of a brane with a vanishing Ricci 
scalar, the limiting perturbations are solely supported by the brane motion in a bulk empty 
from metric perturbations. Moreover, we have also argued that the vDVZ discontinuity 
somehow reappears on a cosmological background where the Ricci scalar does not vanish 
but on a measure zero subset. In this case, if one only considers the time evolution of 
perturbations given by the linearized theory and follows the procedure outlined above for 
fixing initial conditions, one finds divergences in observable quantities at the points where 
the Ricci scalar vanishes when one lets rcH go to infinity. These divergences result in 
a breaking down of the linearized theory, in a similar way to what is happening around 
a static source on a Minkowski background [20,27,30-32], which is intimately related to 
the vDVZ discontinuity. One might be concerned that our results seem to depend on the 
assumption the bulk initial data have a continuous and bounded limit. In fact, the only 
necessary condition is that those bulk initial data are such that the limiting value (as rcH 
goes to infinity) on the brane of the master variable and its derivatives (or of 0, in the 
case of a vanishing are bounded^^, as seen in particular in section l377l When it is 

not the case, we argued that the linear perturbation theory is no longer valid. 
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A Expression of the coefficients C 

Here we give the expressions for the coefficients appearing in equations (|59I63() . 





= l + /x(T + l), 


(1,0) 


= 3{1 + /t(T + 1)}, 
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(2,0) 


= -3/t(1 + T), 
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^(0,1) 
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pp 
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pP 
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^(1,1 

pq 

^A(o,o: 
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'"A(1,0 

pq 
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nq 
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(2 + T) {3 (6iJ^ + 5H^H - + HH^ + 9F^t} 
-/t(2 + T), 

3/t(2 + T), 
2 



3/x 

-rj 



(2 + T) \^AH^H + 6ij2 - 2HH + T (2H^H + H"^ - HH^ } 

37?/t(2 + T)2, 

■2 



rih{2 + Tf 



f2 

•'^ {24^^^^ + 20H'^H + 2HH + T(24iJ^ + IOH'^H - H"^ + HH) + 6T^iJ^}, 



--^(2 + T){T(2 + T)H^ + (6 + 7T + T'^)H^H + (7 + 2T)H'^H'^ + - 
H 

-^{3H\2 + T) + 2H}, 

-^(2 + T){3(2 + T)H^ + BH'^H - H"^ + HH}, 
H 

3/t(2 + T), 
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^{0,1) = (2 + T^){2(2 + T)H^H + (6 + T)ij2 - (2 + T)HH}, 

C(M) = -3r//T(2 + Tf, 
Ca(o,o) = ^{(2 + Tf(-3 + 2T)/76i7+(2 + T)F^ + 2(2 + T)2i/5if 

+(28 + 20T + 3T'^)H^HH - (4 + T)HH^H + (2 + T)i7^((-24 + 5X^)^2 
-(2 + T)lf(3)) _ ii-2^4g ^ -^3^ _ 2T2)ii-3 _ + T)^^^ + (2 + T)HH^^^)^ , 

= -^{3(2T(2 + T)3i7^° + (2 + T)2(18 + 25T + 5T2)i78^-3(-6 + T2)i72^4 

-Tij^ + (2 + T)2(4 + 5T + T^)H'^H + (44 + 52T + 21X2 _^ 2,T^)h^hH 

+(6 + lOT + 3T'^)H^H'^H + (2 + T)HH^H - (2 + T)ii'''((-56 - 56T - lOT^ + T^)H^ 

+(2 + T)//^^)) + H*((42 + 39T + - 3T=^)ii'^ + 2(2 + T)H^ - (2 + T)HH^^'>))} , 

J (6) 

C^(i,o) = -^{2 + T){si2 + T)H' + 5H'H-H' + HH}^^^, 

-3(8 + 6T + T'^)H^HH + (4 + T)HH'^H + (2 + T)i/'^((42 + 19T + T^ji^^ 
+(2 + T)i?(^)) + i?2((54 + 30T + 4T2)ij3 _ 2(2 + T)^^ + (2 + T)HH^^^))^ , 

Cgg^ = ^{(2 + T)2i74 + 3(2 + T)ij2ii--2(3 + T)ii'2 + 2(2 + T)ii"if}^^^, 
Cf4,o) = -3/t(2 + T), 

Cg^ = 7y^{(2 + T)3(6 + T)ii'6iJ+(24 + 6T-T2)ii4 + (2 + T)2(10 + 3T)ii'3_H-^ 

-(12 + 8T + T^)HH'^H - (2 + T)2if4((-20 - 2T + T2)if2 + (2 + T)ij(^)) 
-(2 + T)H^{{-6 + 9T + 4T2)ij3 _ 2(2 + T)^^^ + (2 + T)HH^^^)^ 



(6)' 

C^.i) = (2 + T) { (2 + T)2i/4 + 3(2 + T)i?2ij + (g + T)if2 - (2 + T)//^f } , 

Cf2,i) = 377/t(2 + T)2, 

where iJ^^^ denotes the third derivative of H with respect to the brane cosmic time t, and 
we have defined /x by 

-^(^^ = 6 (^H + 2i72^ being the background scalar curvature of the induced metric on the 
brane. 
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